Pressure— induced amorphization, crystal— crystal transformations and the memory 
glass effect in interacting particles in two dimensions 
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We study a model of interacting particles in two dimensions to address the relation between 
crystal-crystal transformations and pressure-induced amorphization. On increasing pressure at very 
low temperature, our model undergoes a martensitic crystal-crystal transformation. The charac- 
teristics of the resulting poly crystalline structure depend on defect density, compression rate, and 
nucleation and growth barriers. We find two different limiting cases. In one of them the marten- 
site crystals, once nucleated, grow easily perpendicularly to the invariant interface, and the final 
structure contains large crystals of the different martensite variants. Upon decompression almost 
every atom returns to its original position, and the original crystal is fully recovered. In the second 
limiting case, after nucleation the growth of martensite crystals is inhibited by energetic barriers. 
The final morphology in this case is that of a polycrystal with a very small crystal size. This may 
be taken to be amorphous if we have only access (as experimentally may be the case) to the angu- 
larly averaged structure factor. However, this 'X-ray amorphous' material is anisotropic, and this 
shows up upon decompression, when it recovers the original crystalline structure with an orientation 
correlated with the one it had prior to compression. The memory effect of this X-ray amorphous 
material is a natural consequence of the memory effect associated to the underlying martensitic 
transformation. We suggest that this kind of mechanism is present in many of the experimental 
observations of the memory glass effect, in which a crystal with the original orientation is recovered 
from an apparently amorphous sample when pressure is released. 

PACS numbers: 64.70.Kb, 61.50.Ks, 61.43.-j 



I. INTRODUCTION 

Several crystalline materials are known to become 
amorphous under changes of the ambient pressure^S^ 
and this may occur under compression or decompres- 
sion. Among the most thoroughly studied cases, due to 
its technological or geophysical importance, we find H2O, 
SiC>2, GeC>2, AIPO4 and other isostructural compounds. 
The occurrence of pressure induced amorphization (PIA) 
is characteristic of (although not limited to) tetrahedrally 
coordinated materials. 

The simple statement that a material becomes amor- 
phous under pressure covers a broad range of different 
experimental situations and possible theoretical interpre- 
tations. The simplest possibility is a crystalline material 
that transforms (upon change of the external pressure) 
to a completely amorphous structure when a well defined 
transition pressure is reached. We will call this possibility 
true PIA (TPIA). This strict case of PIA is by no means 
the only possible one. A reason for this is that typi- 
cally PIA is claimed to occur (on an experimental basis) 
when the obtained sample does not display sharp X-ray 
diffraction peaks. We will refer to these samples as 'X-ray 
amorphous' material. A relatively large amount of local 
order can still be present in an X-ray amorphous. For in- 
stance, a polycrystal with a sufficiently small crystal size 
is an X-ray amorphous material. An eventual transition 
from a single-crystal sample to this kind of polycrystal 
should be considered, on an experimental basis, to be a 



PIA transition. This is a second case of PIA that we will 
call 'weak' (WPIA). 

Experimentally, WPIA or some variants of it may be 
the rule more than the exception. It is seen in many 
cases that amorphization occurs through a sequence of 
transformations in which one or more intermediate crys- 
talline phases can be identified. For instance, a crystal- 
crystal pressure induced transformation was shown to oc- 
cur in cv— quartz between 21 and 30 Gpa 4,5,6 . The new 
crystalline phases, quartz II and a P1\jc phase, may be 
intermediate states in the path toward complete amor- 
phization, occurring at higher pressures. Yet, the final 
amorphous material may well be only an X-ray amor- 
phous material. 

A common characteristic of both WPIA and TPIA 
is that the starting crystalline sample becomes unstable 
upon pressure. This instability does not typically disap- 
pear if temperature is reduced. This is an indication that 
in most, if not all cases, the mechanisms of PIA can be 
discussed in the limit of T — 0. In this limit, all struc- 
tural transformations that may eventually occur are of a 
mechanical nature, i.e., the displacement of the particles 
during the transformation is strictly guided by the ten- 
dency of the system to minimize its mechanical energy. 
In addition to being a simplifying assumption, the con- 
sideration of vanishingly small temperature is conceptu- 
ally important, as in this case structural transformations 
are necessarily related to mechanical instabilities in the 
system. We notice that although we will present all our 
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analysis at T ~ 0, this is in principle applicable also to 
any temperature below the glass transition temperature 
T g of the fluid phase, where diffusion effects are negli- 
gible. We only need to consider in this case that the 
parameters of the model may depend on temperature. 
Since at T < T g diffusive processes are almost totally 
inhibited, any structural transition (in particular, PIA) 
occurring in this temperature range is of a displacive na- 
ture. This raises the question of the relation between 
PIA and martensitic transformations £ In fact, the dis- 
cussion of the similarities between PIA and martensitic 
transformations is one of the main aims of the present 
paper. 

There is also a long standing controversy about 
whether PIA is of a mechanical or thermodynamical 
naturejSiSiiLiSiAiii^ and this in turn is related to the 
two possible mechanisms of melting at higher temper- 
atures: mechanical or thcrmodynamical. 13 Thermody- 
namical melting implies the growth of the fluid phase 
when it starts to be more stable than the solid phase. 
This growth is heterogeneous and occurs in the presence 
of extended nucleation centers, such as grain boundaries 
or even the surface of the sample. If these defects do 
not exist (something that can be achieved in numerical 
simulations by the use of single-crystalline samples with 
periodic boundary conditions), then the fluid phase is 
not able to nucleate until the ultimate mechanical equi- 
librium conditions are violated. This is known as Born's 
mechanical stability criteria, and the melting mechanism 
in this limit is called mechanical melting. Mechanical 
melting is homogeneous, as it affects the whole sample 
at the same time, and occurs very rapidly once the sta- 
bility limit is reached. It has been suggested that these 
two melting mechanisms can be extended to describe two 
different amorphization mechanisms. 13 

There are quite a few numerical models devised to 
study PIA in different materials. In some cases, a thor- 
ough description of a particular material is pursued. In 
this case, the models range from first principle density 
functional calculations^^ to molecular dynamics simu- 
lations with different pair-potential interactions>iSiiLi& 
Within these models, PIA can be directly studied in 
concrete systems such as H2O or SiC>2. A more 
qualitative approach is obtained by using core-softened 
potentialSfi2i2£ which proved to be very useful for the 
study of water-like anomalies and polyamorphism in 
tetrahedrally coordinated materials 2L2i/S^iS^i2k The aim 
of this approach is not to carry out an analysis of the 
properties of a specific material, but to find out general 
trends and qualitative behaviors that shed light on the 
properties of a large family of materials. This is the ap- 
proach that we will follow here. 

We study here a two-dimensional model of identical 
particles interacting through a spherically symmetric two 
body potential. We have recently shown that choos- 
ing appropriately the interaction potential, the model 
presents TPIA upon pressure release if we start from the 
stable high pressure phased Here we analyze in detail 



the WPIA that occurs when we increase pressure starting 
from the stable low pressure phase. We show that this 
transformation is of a martensitic nature, and analyze 
the nucleation and growth of the crystallites of the new 
phase. Depending on details of the potential, two differ- 
ent extreme possibilities occur. In one of them, the new 
phase grows in large crystallites of the (three) different 
martensitic variants, and crystallites grow perpendicu- 
larly to the habit line of the martensitic transformation. 
This is a case of ideal martensitic transformation. In the 
other limiting case, the new crystallites nucleate but can 
hardly grow. In this WPIA scenario is realized. 

The final sample is X-ray amorphous, but it conserves a 
memory of its original starting configuration. This mem- 
ory manifests upon pressure release: the sample recovers 
(partially) the orientation it had prior to compression. 
We discuss the characteristics of the interaction poten- 
tial that are responsible for the different behavior. 

The paper is organized as follows. In Sec. [n] we 
present the model and details of the performed numer- 
ical simulations. The ideal martensitic transformation 
is introduced in Sec. IIIII In Sec. IIVI we describe the 
martensitic transformation for the case when low energy 
barriers are involved in the growth process. The case 
were this energy barriers are higher, which leads to a 
WPIA, is discussed in Sec. Finally, Sec. IVII contains 
some discussion and conclusions. 



II. THE MODEL, NUMERICAL DETAILS AND 
EQUILIBRIUM PROPERTIES 

We consider a two-dimensional system of particles, in- 
teracting through a specially devised isotropic, purely 
repulsive pairwise potential. Its main characteristic 
is the existence of two competing distances tq and 
7"i > tq at which neighbor particles prefer to be lo- 
cated, depending on the applied pressure. This kind 
of core-softened potential has been previously used to 
study anomalous properties of tetrahedrally coordinated 
materialsjSLSSiSi^ and it was recently used to study the 
TPIA case also.— The interaction potential V(r) between 
two particles separated a distance r is explicitly given 

J-. ..24.26 

V(r) = -^-[8(r/a-1.55) 2 + 16(r/a-1.55) 4 + 
0.81 

w{r/a- 2) 2 + 0.81] for r < 2a (1) 
= for r > 2a 

Here, e and a fix the energy and length scales, and w is 
a control parameter that allows us to change the height 
of the shoulder of the interparticle potential (see Fig. 
Neighbor particles prefer to be located on both sides of 
the shoulder, in the positions qualitatively indicated in 
Fig. \T]as r o an d r i- Simulation results will be presented 
mainly for w — 6 and w = 10 to assess the effect of the 
barrier height on the morphology of the microstructure 
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obtained. Note that, although for w = 10 the curva- 
ture of the interparticle potential does not change sign, 
the potential gives rise to two preferred distances when 
analyzing two- or three-dimensional configurations. 

The system is simulated by standard molecular dy- 
namics in the NVT ensemble i2L22i This choice allows us 
to survey all regions of the volume-pressure curve, in- 
cluding those that would be unstable in constant pressure 
simulations. We perform simulations with 9540 particles, 
with a time step of St — O.Olio (where to = a^fmfk is the 
time unit, m being the mass of each particle). Pressure 
and elastic constant are calculated along the simulations 
using standard formulas (see the Appendix for details). 
We always use periodic boundary conditions (PBC). This 
is a particularly convenient choice when we want to avoid 
the heterogeneous nucleation of the new phase. On the 
contrary, in the case in which we study the heterogeneous 
nucleation, we introduce the nucleation centers by hand 
in the crystalline structure. The use of PBC in this case 
does not impose any strong additional restriction in the 
possible transformations to be observed, as our systems 
contain enough atoms to accommodate essentially any 
possible new phase, even using PBC. 

The volume of the system is changed at a fixed rate 
during the simulation. This change is done by rescaling 
all coordinates of the particles and size of the simulation 
box. The volume change rate is taken to be as low as 
possible (within reasonable computational time) to sim- 
ulate a quasistatic process. However, it is of course a 
very large rate in physical units. For instance, a typical 
volume change rate we use is lO" 4 ^ 1 , and for an atomic 
particle (to ~ 10 _15 s) this implies a rate of the order of 
10 -2 per femtosecond. 

By rescaling down the velocities of the particles dur- 
ing the simulation, temperature is kept constant at a very 
small value of the order of 10 _5 e (we take the Boltzmann 
constant to be unity) . As the transformation process im- 
plies local mechanical instabilities, there is a systematic 
transformation of potential energy into kinetic energy. 
In experiments, this reflects an increase of the temper- 
ature of the system. Here the kinetic energy excess is 
taken away through the mentioned velocity rescaling pro- 
cedure. 

We emphasize that the use of a two-dimensional system 
is dictated by the simplicity of the simulation and also 
of the interpretation of the results. While we expect the 
same qualitative kind of results to reappear in full three 
dimensional systems, the computational time increases 
noticeably, and the adequate identification of the phases 
becomes much more involved. Note that in a previous 
work it has been shown that the number of equilibrium 
solid phases noticeably increases when passing from two 
to three dimensions^ then adding to the complexity of 
the three-dimensional case. 

One of the most striking properties of particles inter- 
acting through the potential described in the previous 
section is that the stable crystalline ground state is not 
necessarily the triangular lattice. In fact, although at low 




FIG. 1: Pairwise core-softened interparticle potential used in 
the present work [from Eq. Q]. The distances ro and n are 
schematically shown. 

(high) pressures the ground state of the system is a trian- 
gular structure with a lattice parameter around r\ (ro), 
at intermediate pressures lower energy configurations can 
be obtained that have some of its first neighbors at dis- 
tances ~ ro , and some others at ~ r\ . By increasing the 
pressure, the configuration to which the low density tri- 
angular (LDT) structure is expected to transform is the 
chainlike (CH) structured This is the transition that we 
will be mostly interested in, and we will analyze it for 
two different values of w, namely w = 10 and w = 6. We 
will see that there are noticeable differences between the 
two cases. 



III. THE IDEAL MARTENSITIC 
TRANSFORMATION 

We will study the behavior of the LDT phase upon 
volume reduction. The LDT-CH transformation admits 
a natural description as a martensitic transformation as 
depicted in Fig. [21 In this description, the LDT corre- 
sponds to the austenite phase, and the CH to the marten- 
site phase. Upon volume reduction, the martensite phase 
becomes stable above the equilibrium transition pressure 
P e . Two structures can coexist at this pressure, namely 
a LDT structure with lattice parameter r e [Fig. Ufa)], 
and a CH structure (corresponding to an oblique two- 
dimensional Bravais lattice) with parameters r$, rf, and 
j3 e [Fig. Hfb)]. To minimize the interface energy at the 
equilibrium pressure, the two structures have to match 
along a particular line, called the invariant line or habit 
line (habit plane in three dimensions, line AB in Fig. 
121 . This line is determined from the condition that both 
phases are unstrained away from iti A relative rotation 
of the two phases is necessary in order to define the in- 
variant line. The orientation of the invariant line can be 
easily evaluated given the parameters of the martensite 
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TABLE I: Structural parameters at the LDT-CH equilib- 
rium point (see Fig. I^Jfor definitions of the crystallographic 
parameters) , along with the pressure and volume values at the 
mechanical instability. The angles 9 e and 4> e correspond to a 
relative rotation of the two phases and a reference angle be- 
tween the invariant line and the LDT structure, respectively 
(see Fig. gj. 



Parameter 


w = 6 


w = 10 


P e a A /e 


4.83 


8.55 


v e /a 2 (LDT) 


3.01 


2.74 


r e /a(LDT) 
^ e /a 2 (CH) 


1.86 


1.78 


2.29 


2.38 


r 7«(CH) 


1.25 


1.37 


rf/a(CH) 


1.93 


1.87 


/3 e (CH) 


71.07 


68.56 


e (CH) 


1.46 


1.65 


(CH) 
P c a 2 /e 


5.47 


9.43 


6.38 


9.34 


v c /a 2 


2.76 


2.64 



and the austenite phase. The structural parameters of 
the austenite (LDT) and martensite (CH) phases at the 
equilibrium pressure are given in Table [I] 

The ideal transformation between the LDT and CH 
structures does not imply bond breaking, which means 
that each particle conserves the same first neighbors in 
the two phases. This implies microscopic memory in the 
sense that each local neighborhood transforms back to 
exactly the original configuration upon decompression. 
Furthermore, this implies that after a cyclic LDT-CH- 
LDT transformation, the original structure is recovered 
with the same initial orientation. This property is di- 
rectly related to the shape memory effect of martensites. 

If an invariant line separating the two coexisting phases 
has formed at the transformation pressure, the transfor- 
mation can proceed as volume is reduced by reaccom- 
modating particles of the austenite close to the inter- 
face to fit into the sites of the martensite^ The advance 
of the interface occurs through the lateral movement of 
the steps that form the interface, and at the equilibrium 
pressure it implies the surmounting of an energy bar- 
rier. Then, at T — a small overpressure is necessary 
to climb this energy barrier. As we will see, this ideal 
growth mechanism in the direction perpendicular to the 
invariant line becomes more complicated when strains 
are involved, and depends strongly on the height of the 
barrier for the transformation between the austenite and 
martensite phases. 



IV. NUCLEATION AND GROWTH IN THE 
CASE OF A LOW— ENERGY BARRIER 

The previous ideal description of the martensitic trans- 
formation can differ from the actual way in which the 
martensite nucleates and grows from the pure austenite 
phase in a real situation. In this section, we study this 




FIG. 2: Ideal (a) LDT and (b) CH structures. The lat- 
tice parameter and the crystallographic directions are shown, 
together with the perfect and distorted hexagonal unit, (c) 
Interphase between the LDT and CH structures at equilib- 
rium. The invariant line AB is rotated (j> e with respect to 
the LDT structure. The [1,0] direction of the CH structure 
is rotated 9 e with respect to the [1,0] direction of the LDT 
structure. 



process using the potential in Fig. ^with w = 10- On 
increasing pressure, the low density triangular structure 
becomes thermodynamically unstable against the chain- 
like structure beyond P e . However, as this transforma- 
tion is first order, at T — the perfect LDT structure 
survives in a metastable state at higher pressures until 
it finally becomes mechanically unstable. For the case of 
a perfect (without defects) starting sample, we call the 
instability pressure P c (and the corresponding specific 
volume v c ). The mechanical stability limit corresponds 
to the value of pressure at which one of the normal modes 
(phonons) of the original structure first destabilizes. The 
calculation of the phonon structure of the triangular lat- 
tice is straightforward. For our model potential Eq. 
interactions with second and further neighbors are iden- 
tically zero, and the stability condition can be written 



d 2 V 3 dV 



dr 2 



dr 



>0, 



(2) 



where r is the nearest neighbor distance of the triangular 
lattice. At the instability pressure P c , three equivalent 
shear phonon branches become zero energy. They corre- 
spond to shear phonons in the three equivalent directions 
[0,1], [1,0], and [-1,1] (Fig. [2$. Note that, due to the 
lack of next nearest neighbors interactions, the phonon 
branch along these directions is of the type ~ sin(fc), and 
at P c the whole phonon branch becomes dispersionless. 

The fact that the destabilization is initiated by a shear 
phonon is consistent with an analysis in terms of Born 
inequalities 31 . In fact, condition © is equivalent to 
the generalization of Born stability criteria to stress de- 
pendent conditions 32 ' 33 in the static (zero temperature) 
limit. In two dimensions, the Born stability criteria es- 
tablish that a system is stable if the stress dependent 
shear modulus G' — (Cn — Ci 2 )/2 — P is positive^ 3 -. For 
the zero temperature case G' = y/3/4:[V"(r) + 3V'(r)/r] 
(see Appendix), and then the Born stability criteria re- 
duce to condition 

The results of numerical simulation [Fig. GJa)] con- 
firm the previous analysis. Starting from a perfect tri- 
angular structure at zero temperature, once the system 
reaches the instability volume v c , the triangular configu- 
ration destabilizes and a new stable configuration should 
emerge. The evolution of the stress dependent shear 
modulus upon volume decrease is shown in Fig. l3TK). 
where it is clear that it vanishes at the mechanical insta- 
bility point. 

The consideration of a perfect sample without defects 
is important in order to realize the equivalence of the 
mechanical instability limit with the vanishing of the en- 
ergy of a phonon mode. Moreover, in the present case it 
also coincides with the vanishing of an elastic constant (a 
long wave length elastic constant), because of the lack of 
interactions to next nearest neighbors^ But certainly, a 
defect free sample is not typically encountered in experi- 
ments. It is important to realize that even tiny amounts 
of defects can produce important changes in the trigger- 
ing of the transformation. For instance, a single vacancy 
can destabilize the lattice at pressures much lower than 
P c . This was shown already in the case of TPIA^& and 
it is also true here. The evolution of pressure upon vol- 
ume decrease for a lattice with a single vacancy is also 
shown in Fig. EUa). We see that a single vacancy is able 
to destabilize the original triangular lattice without es- 
sentially any pressure overshooting. Note, however, that 
the inclusion of a single vacancy has no appreciable effect 
on the values of the shear modulus prior to the transition 
[Fig. 13(b)] . This means that even in the presence of tiny 
amounts of disorder, the examination of the clastic con- 
stants (or even of the whole phonon spectrum) will not 
indicate the approach to the instability point. In fact, 
the instability occurs because of the destabilization of 
localized modes close to the defect, which have no effect 
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FIG. 3: (a) Volume- pressure curves for w — 10. The dotted 
line is the analytic expression for the evolution of the LDT 
structure if this remains stable, and the thin continuous line 
marks the equilibrium transition to the CH structure and the 
subsequent evolution of this structure. The full and open 
symbols correspond to simulations with and without a va- 
cancy in the original sample. Vertical dashed arrows indicate 
the compression-decompression paths, (b) Evolution of the 
stress dependent shear modulus G' = (Cn — (7i2)/2 — P upon 
decompression for w = 10. The dotted line is the analytical 
result for G' evaluated for the perfect triangular lattice (see 
the Appendix). Note how in the defect free case, the approach 
to the transition point (v —* v c ) is signaled by the vanishing 
of the stress dependent shear modulus. 



on the spectrum of the (extended) phonons. The absence 
of macroscopic effects signaling the approach to the in- 
stability is usually associated to thermally activated first 
order transitions. We see that a similar phenomenology 
occurs here in the absence of thermal fluctuations, and 
within a scenario of a transition driven only by mechan- 
ical instabilities. 

In Fig. a sequence of snapshots is shown to follow 
the morphology of the transformed regions in the case 
in which a single vacancy is present in a monocrystalline 
sample. We can see that the first instability produces the 
collapse of adjacent rows of particles, around the defect. 
These collapses generate three platelets that propagate 
out of the vacancy [Fig. Ufa)]. The collapse directions 
in the three platelets are different, and are related to the 
three different variants of the martensitic transformation. 
In this stage, the martensite crystallites are compressed 
in the longitudinal direction, in order to match the lat- 
tice parameter of the original triangular structured In 
a subsequent stage [Fig. 0(b), (c)], the new crystallites 
reaccommodate to take the form of wedge shaped crys- 
tallites, with the interface between the martensite and 
austenite phase being very close to the ideal habit line of 
the transformation. Then upon further compression [Fig. 
Uc-e)] the crystallites grow perpendicularly to the habit 
line, very much as described for the ideal case in the pre- 
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FIG. 4: (Color online) Snapshots of the system with w — 10 
corresponding to: (a) v/a 2 = 2.7, (b) v/a 2 — 2.6, (c) v/a 2 — 
2.5, (d) v/a 2 = 2.4, (e) v/a 2 — 2.3, on compression, and (f) 
v/a 2 — 2.8 after decompression from the configuration in (e). 
Only collapsed particles are shown joined by segments, with 
three different colors to identify the three different variants 
of the martensite (intact austenite is not plotted for clarity) . 
The starting system contains a single vacancy at the position 
where the three different variants start nucleating in (a) . Note 
that periodic boundary conditions are used. 



vious Section. The final state of the system (when the 
whole sample has transformed to the martensitic phase) 
consists of a polycrystalline sample with the three differ- 
ent orientation of the possible variants [Fig. [He)]. Note 
from Fig. Eta), that after nucleation, pressure has re- 
mained essentially constant during the whole transfor- 
mation from the LDT to the CH structure, at a value 
slightly larger than the equilibrium value. This value re- 
flects the overpressure needed to keep the transformation 
proceeding. 

If the volume of the system is increased back start- 
ing from the sample completely in the martensitic phase, 
upon complete pressure release essentially any particle 
in the system returns to its original position, in which it 
has exactly the same neighbors as before compression. In 
particular, this implies that the original monocrystalline 




FIG. 5: Evolution of the averaged structure factor S(k) with 
decreasing volume for the LDT-CH transformation with w = 
10 (the curves are vertically displaced, for clarity). Different 
curves correspond to the configurations in Fig. El( a )~( e )- 





FIG. 6: Two-dimensional (upper) and angularly averaged 
(lower) structure factors of the configuration in Fig. EJe), and 
of the configuration obtained by quenching from the melt at 
the same density. In the top, the two-dimensional structure 
factors of the polycrystalline CH structure (left) and quenched 
amorphous structure (right) are shown. In the lower figure, 
bold and thin lines correspond to the averaged S(k) for the 
polycrystalline CH and quenched amorphous structures, re- 
spectively. There is a clear crystalinity observable in the first 
case, even if we look only at the radial structure function. 
This sample should be characterized as polycrystalline. 
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austenite sample is reobtained, with the same orientation 
as in the starting configuration^ This is not a surprise, 
since it is in the essence of the shape memory effect of 
martensites. It can be seen in Fig. Ola) that the re- 
turn occurs with a noticeable hysteresis with respect to 
compression. This hysteresis is a consequence of the indi- 
vidual hysteresis of pairs of neighbor particles on passing 
between the two equilibrium distances, and can be used 
to quantify the energy dissipated during a compression- 
decompression cycle. 

Now we want to discuss how the transformation man- 
ifests when we observe the diffraction patterns of the 
structures. In Fig. OD we show the angularly averaged 
structure factor S(k) corresponding to the configurations 
in Fig. EJa)-(e). It is possible to see how some character- 
istic peaks of the LDT structure progressively disappear 
while new peaks appear that correspond to the CH struc- 
ture. A broadening of the peaks is also apparent, due to 
the finite size of the crystallites in the new configura- 
tion. This broadening however is not large enough to 
consider the sample as X-ray amorphous. By comparing 
the structure factor of this structure with that of a sample 
obtained by quenching from the melt at the same density 
(Fig. |SJ, we first observe that the two-dimensional struc- 
ture factor displays clear evidence of the polycrystallinc 
character of the sample. But even if we look only at the 
radial structure factor, there are much sharper peaks for 
our sample compared to the quenched liquid. Then the 
configuration in Fig. ^e) can be correctly identified from 
its structure factor as a polycrystalline sample. 



V. BLOCKED GROWTH FOR LARGE ENERGY 
BARRIER 

In the present section, we will study the modification 
brought upon by a change of the parameter w of the po- 
tential. The value used here is w = 6. In principle, the 
same description of the ideal martensitic transformation 
is valid also in this case. The value of the crystallographic 
parameters are given in Table [I] In Fig. EI a), we show 
the volume-pressure evolution of the volume controlled 
compression for w = 6. Open symbols correspond to the 
starting sample without defects. The simulation shows 
that the system remains in the triangular structure un- 
til the mechanical stability limit (P c ,v c ) is reached. On 
the contrary, the sample with a single vacancy begins to 
transform when local mechanical instabilities occur, be- 
fore reaching the mechanical stability limit. The behav- 
ior of the stress dependent shear modulus for the present 
case, with and without the vacancy, is shown in Fig. I3b). 
The phenomenology up to here is very similar to the case 
of the previous section. However, important differences 
appear when we look at the morphology of the samples 
and the structure factors. We again focus on the case 
of a system with a single vacancy. In Fig. [8] snapshots 
of the system as volume is decreased are shown. As in 
the previous case, the first elements of the new phase to 
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FIG. 7: (a) Volume-pressure curves for w — 6. The dotted 
line is the expected evolution of the LDT structure, and the 
thin continuous line marks the equilibrium transformation to 
the CH structure and the subsequent evolution of the new 
structure. The full and open symbols correspond to simula- 
tions with and without a vacancy in the original sample. Ver- 
tical dashed arrows indicate the compression-decompression 
paths, (b) Evolution of the stress dependent shear modulus 
G' = (Cii — Cvi)j2 — P upon decompression for w = 6. The 
dotted line corresponds to G' evaluated from the interatomic 
potential (see the Appendix). P c and v c are given by the 
vanishing of the stress dependent shear modulus in (b) (see 
TableQ. 



be nucleated are the platelets starting at the position of 
the vacancy. But contrary to the previous case, they do 
not rearrange to generate an interface along a habit line. 
As the growth process in these conditions implies much 
higher overpressure, the growth remains mostly blocked, 
and upon further compression nucleation of new platelets 
throughout the sample occurs. New and old platelets in- 
teract elastically, and deform each other, generating at 
the end a structure that in spite of being martensitically 
related to the original one, contains crystallites that are 
extremely small. 

On very general grounds, since our polycrystalline 
structure has very small crystallites, a broadening of the 
characteristic peaks of the azimuthally averaged X-ray 
pattern is expected^ There is a continuous crossover 
however from this 'broadened peaks' scenario, to one in 
which no sharp peaks are observed. In the last case the 
sample should be characterized as X-ray amorphous. But 
even in this case it can possess specific crystal character- 
istics, and in particular it may conserve a structural re- 
lation with the parent structure that will manifest upon 
decompression. 

In Fig. we show angularly averaged structure fac- 
tors corresponding to the configurations shown in Fig. |S| 
(a)-(e). Again, the peaks of the LDT structure disappear 
while reducing the volume, but the intensity of the peaks 
corresponding to the CH structure is much weaker than 
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FIG. 8: (Color online) Snapshots of the volume controlled 
simulation of the system with to = 6 corresponding to: (a) 
v/a 2 = 2.95, (b) v/a 2 = 2.8, (c) v/a 2 = 2.6, (d) v/a 2 = 2.3, 



(e) v/a 2 



2.1, 



on compression, 

.2 



and (i) v/a = 3.05 on 
decompression from u/a 2 = 2.1. Only collapsed particles 
are shown joined by segments, with three different colors to 
identify the three different variants of the martensite (intact 
austenite is not plotted for clarity). The starting system con- 
tains a single vacancy at the position where the three different 
variants start nucleating in (a) . Note that periodic boundary 
conditions are used. 



in the w = 10 case. In Fig. ^]the two-dimensional X- 
ray diffraction patterns of the structure obtained from a 
single crystal LDT sample, and an amorphous structure 
obtained by quenching from the melt (with the same den- 
sity), are shown. We see that the two-dimensional pat- 
tern of the sample obtained by compression possesses an 
angular structure reminiscent of the original crystal ori- 
entation (in fact, this structure will be responsible for 
a 'memory effect' when the sample is decompressed, see 
below). However, if our starting sample is already poly- 
crystalline, as is usually the case experimentally, the X- 
ray pattern is by definition isotropic. In that case all the 
information available is the radial X-ray pattern. The an- 
gularly averaged S(k) are shown at the bottom of Fig. ITU1 
On the basis of the S(k) function, the compressed sam- 
ple cannot be clearly distinguished from that obtained 




FIG. 9: Evolution of the averaged structure factor S(k) with 
decreasing volume for the LDT-CH transformation with w = 
6. Each curves correspond to the configurations in Fig. |H1 
(a)-(e). 



by quenching from the melt. This sample is then X-ray 
amorphous. Thus, we have here a sample that has WPIA. 

The fact that this sample is X-ray amorphous, but still 
has within its structure some information of the parent 
crystalline phase from which it was obtained may be re- 
lated to the so-called memory glass effect (MGE). The 
MGE was first identified in AIPO4 berlinite) 37 ! 38 and 
then it was recognized to occur in a variety of materi- 
als such as SnLj and LiKSO^ The effect consists in the 
amorphization of a sample under pressure, and the suc- 
cessive recovery of a crystalline sample, with the same 
original orientation upon decompression. In the follow- 
ing, some characteristics of the MGE for the present 
model system are studied. 

The complete v — P curves during volume controlled 
compression-decompression cycles are shown in Fig. 1111 
for different values of the final volume reached, v m i n . 
The final configurations after the cycling are shown in 
Fig. All the snapshots correspond to a recovered 

volume of v/a 2 = 3.05. Shown are the defects intro- 
duced by the cycling, corresponding to small regions of 
the different martensitic variants. Small dots correspond 
to the position of the particles. On a direct examination 
of the recovered structures, it is seen that the system has 
a tendency to recuperate its original orientation upon 
decompression. A closer inspection shows that the re- 
covered LDT structure retains its initial orientation for 
high v m i n , but becomes polycrystalline for lower v m i n , 
losing the 'memory' of the initial orientation. 

It is convenient to have a single parameter that mea- 
sures the degree of memory of the structure. We take 
this parameter to be W, defined as 



W 



36N 



k,k' 



(3) 



9 




k a 



FIG. 10: Comparison of the structure factor corresponding 
to the sample obtained in the LDT-CH transformation [Fig. 
ISfe)] with a quenched amorphous structure at the same den- 
sity. In the top the two-dimensional structure factors for the 
polycrystalline CH structure (left) and quenched amorphous 
structure (right) are shown. The angular averaged S(k) is 
plotted at the bottom for the polycrystalline CH structure 
(bold continuous line) and quenched amorphous structure 
(thin continuous line). 
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FIG. 11: v — P compression-decompression cycles for differ- 
ent values of the final volume reached. From bottom to top 
Vmin/a 2 goes from 1.6 to 2.6, in steps of 0.1. Dotted and thin 
continuous lines as in Fig. [7| 



where k and k' are the nearest neighbors of particle j be- 
fore and after the compression-decompression cycle, and 
ajk(k') is the- angle between ?^W) ~ an d an arbitrary 
reference direction. In particular, it can be seen that if 
the sample returns exactly to its initial configuration, we 
get W = 1. On the other hand, if the final local orien- 
tations in the sample are completely uncorrelated with 
the original ones, then W = 0. We take as the reference 
volume to compute W the value v/a 2 — 3.05 and the plot 
of W as a function of the minimum volume v m i n reached 
during compression is seen in Fig. 1131 We see in fact, 
that this number is close to one for v m i n /a 2 > 2, whereas 
it is essentially zero when v m i n /a 2 < 1.9. The memory of 
the system is lost continuously when v m ,„ becomes lower, 
but note the interesting fact that for v m i n =2.1, in which 
there is no LDT phase left in the compressed sample (see 
Fig. [5J , the value of W is still clearly different from zero. 
All the memory of the system here can be associated to 
the preferred orientation of the crystals of the CH phase, 
which in turn occurs because of the martensitic nature of 
the LDT-CH transition. As for this volume the sample 
is X-ray amorphous, this is an example of MGE. 



VI. DISCUSSION AND CONCLUSIONS 

We have studied a two dimensional model at T = 
that displays a pressure driven martensitic transforma- 
tion between a triangular structure (the austenite) with 
rotational Cq symmetry, and a chain like structure (the 
martensite) with a lower (C2) rotational symmetry. The 
martensite phase can be considered to be obtained from 
the austenite through a compression in one direction and 
an expansion in the perpendicular direction. The trans- 
formation occurs with an appreciable volume change. 
This is natural since pressure (which couples to volume in 
the free energy) is the driving force of the transformation. 

The transition does not spontaneously occur at the 
equilibrium transformation pressure because an energy 
barrier exists. An excess pressure is needed to trigger the 
transformation. Defects in the original crystalline struc- 
ture produce the appearance of 'reaction paths' which 
have (at the equilibrium pressure) relatively low energy 
barriers. These barriers vanish with a relatively low 
overpressure, triggering the transition. The necessary 
overpressure in the presence of defects is typically much 
smaller than that needed to destabilize the defect free 
structure. We have shown in particular, that vacan- 
cies are very effective in favoring the nucleation of three 
platelets related to the three variants of the martensite, 
which are however elastically distorted with respect to 
the ideal martensite. We call them pseudo-martensitic 
crystals^ 

After the nucleation of pseudo-martensite platelets, 
two different evolutions were observed. In one case, 
the pseudo-martensitic platelets find the way to become 
wedge shaped crystallites of the true martensitic struc- 
ture. In this process, the internal stress of the martensite 
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FIG. 12: (Color online) Snapshots of the system as recovered from the compression-decompression cycles shown in Fig. 1111 All 
the snapshots correspond to a final volume of v/a 2 — 3.05. The minimum volumes reached are (a) v/a 2 — 2.4, (b) v/a 2 = 2.2, 
(c) v/a 2 = 2.0 and (d) v/a 2 — 1.8. Different colors indicate particles collapsed along each of the directions of the possible 
martensitic variants. Dots corresponding to the position of the particles are shown in order to identify the appearance of 
polycrystallinity. 



relaxes, and the interface between austenite and marten- 
site becomes very nearly the ideal invariant line of the 
transformation. These wedge shaped crystals grow upon 
volume decrease, until the whole sample consists of a col- 
lection of crystallites of the three different variants of the 
martensite, and no austenite remains. In this case, the 



final size of the crystallites is expected to depend mainly 
on the density of defects in the system and the compres- 
sion rate, since there are no other important limitations 
to martensite growth. In principle, very large crystal 
sizes can be obtained in samples with low concentration 
of defects driven at very a low compression rates. 
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FIG. 13: The parameter W defined in Eq. © as a function 
of the minimum volume v m in reached in the compression- 
decompression cycles shown in Fig. 1111 



In a second case, the conversion of pseudo-martcnsitc 
to true martensite is strongly inhibited. In this case, 
upon compression, new platelets of pseudo-martensitic 
structure nucleate, that interact with preexisting ones. 
The structure obtained when the austenite has com- 
pletely disappeared is much more disordered than in the 
previous case, and in particular, the final crystallite size 
in this case is expected to remain of the order of few 
interatomic distances. In some cases, the angularly av- 
eraged diffraction pattern of the final structure may well 
be taken by that of an amorphous structure. We have 
called this phenomenon weak pressure induced amor- 
phization. This is not a true pressure induced amorphiza- 
tion since the final structure has a 'memory' of its parent 
crystal that is observable, for instance, in the full two- 
dimensional diffraction pattern. Upon decompression of 
this sample, a defective crystal is obtained, which, how- 
ever, has a consistent tendency to be oriented in the same 
form as the original parent crystal. This is a description 
of the 'memory glass effect' (MGE), observed in different 
experimental situations. 

The two different behaviors were obtained by changing 
the parameters of the interparticle potential. Concretely, 
the easy growth of the martensite was observed when the 
barrier between the two possible interparticle distances 
was small, whereas inhibited growth was observed when 
this barrier was higher. 

Although the present model is not aiming to repro- 
duce the characteristics of a specific real system with its 
complex interactions, some similarities concerning PIA, 
martensitic transformation and the MGE, may be ad- 
dressed between this simple model system and real sys- 
tems. 

As we neglect temperature effects in our molecular 
dynamics simulation, the transformation between the 
austenite and martensite phases is of a mechanical nature 



(and this should be qualitatively correct for any T < T g ). 
This means that the austenite is driven into a metastable 
region, until a mechanical instability limit is reached, and 
the martensite starts to nucleate spontaneously. Note the 
important fact (already present in our previous results 
inS 6 .) that once nucleated, the interface between the two 
phases can propagate spontaneously, and this propaga- 
tion, which occurs essentially at the sound velocity, is not 
limited by the diffusivity of the particles in any of the two 
bulk phases (which is zero here). In the case of defect free 
samples the mechanical instability of the austenite is sig- 
naled by the softening of a phonon mode. However, in the 
presence of defects the first mechanical instability corre- 
sponds to localized modes around the defects, and there- 
fore it is not observable in the phonon spectrum. Thus, 
our present and previous 2 ^ results reinforce our view that 
PIA is a mechanical transition, i.e. always triggered by 
mechanical instabilities. We find that defects can cause 
the transition pressure to be depressed to a value very 
close to the ideal thermodynamical equilibrium pressure 
between the two phases. We recognize however, that this 
result can be not universal but system dependent. The 
relevance of defects density in lowering the amorphiza- 
tion pressure was recently experimentally highlighted^ 
although whether the amorphization pressure approaches 
the extrapolated melting line upon introduction of more 
defects is not clear. The experimental evidence of amor- 
phization occurring near the extrapolated melting line in 
some systems is at the base of a 'thermodynamical' sce- 
nario for PIA as directly related to the 'thermodynamical 
melting' mechanism at higher temperatures. We see that 
this evidence is also fully compatible with a mechanical 
mechanism. 

A transformation reminiscent of our martensitic trans- 
formation with easy growth occurs between the high pres- 
sure phases of GaPOzi. In this case, the low-cristobalite 
phase of GaP04 transforms to a Cmcm structure on pres- 
sure loading. 4 Sdl The transformation is of a displacive 
nature, and the mechanism consists of a fourfold-to- 
sixfold coordination change of the gallium coordination 
shells The peaks of the polycrystal are clearly identi- 
fied, which indicates a large crystal size. In view of the 
results presented here, it could be interesting to exploit 
the martensitic point of view in analyzing this kind of 
pressure induced transformation. 

On the other hand, the transformation in the case 
of blocked growth is in line with investigations of 
the high pressure transformations of AIPO4 and the 
MGEi 42 i 43 i 44 i 45 i 46 What was first viewed as a true 
PIA with MGE, was later identified as a crystal- 
crystal pressure-induced transformation between a berli- 
nite AIPO4 structure and a Cmcm polycrystalline struc- 
ture with very small crystal size, i.e., in our language, a 
WPIA. Carefully done molecular dynamics simulations 43 
and X-ray experiments^ have identified some peaks of 
the X-ray diffraction patterns of the high pressure phase 
as corresponding to the Cmcm structure. In addition, re- 
cent work has been devoted to establish that the Cmcm 
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structure is the stable one at high pressures^ but lit- 
tle is known about the mechanism of the transformation. 
In view of the similarities between the MGE in our two- 
dimensional model and in the AIPO4 case, we suggest 
that MGE occurs simply as a consequence of the amor- 
phization being only of the weak type, instead of a true 
one. Furthermore, we speculate that the mechanism of 
the berlinite- Cmcm transformation in AIPO4 could be 
very similar to the mechanism of the low cristobalite- 
Cmcm transformation in GaPC>4. However, the growth 
of the Cmcm structure over the berlinite phase of AIPO4 
could be inhibited by energetic barriers, leading to the 
polycrystal with very small crystal size. 

To conclude, we notice how the present model, al- 
though highly simplified and not representing any par- 
ticular real system, captures much of the phenomenology 
observed in a variety of real cases. As such, it is an im- 
portant theoretical tool to guide the analysis of concrete 
cases. We expect to be able to go further in this direction 
in the future. 
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APPENDIX: STRESS DEPENDENT ELASTIC 
CONSTANTS AND THEIR STATIC LIMIT 



x r^apr^apr^apr^ap) . (A. 3) 

Here, V'(r) and V"(r) correspond to first and second 
derivatives of the interaction potential energy at r. 

Consider the static (zero temperature) limit. In this 
limit one has to consider only the perfect crystal struc- 
ture without any thermal fluctuation. In this case the 
first term of Eq. (| A. 1|) . proportional to the temperature, 
vanishes. The second term, which measures the fluctu- 
ation of the pressure tensor, vanishes too. Besides, the 
term containing the momentum of the particles in Eq. 
(IA.2II is zero, and Cijki — B>ijki- Then, in this limit, the 
calculation of the elastic constant is straightforward, and 
instead of using the complete fluctuation formulas, which 
needs an ensemble average to be performed, we use the 
static limit Eq. (| A. 3|> calculated over a given sample. 

Let us consider, for instance, the calculation of C\i in 
a triangular structure of lattice parameter r with nearest 
neighbors interactions only. 




Fluctuation formulas for stress-free elastic constants in 
the NVT ensemble are given byS 



2NkT 



ijkl 



V 

({Pij Pu > - (Pij ) (Pa ))+B m , ( A. 1 ) 



where N is the number of particles, k is the Boltzmann 
constant, T is the temperature and V is the volume of the 
system. The brackets, (...), denotes ensemble averages. 
Pij are the components of the pressure tensor given by 



P — 



P(i) a P(j) a v-^ V'(r aP ) 



E 

a</3 



r a /3 



(A.2) 

with a and (3 particle running indexes, P(j,) a the momen- 
tum of the particle a (with mass m a ) in the i direction, 
and r a @ the distance between particles, with component 
r (i)af3 m the * direction. Finally, 



a<0 \ 



V"{r a p) V'{r af} ) 



'a/3 



where we used N/V = 2/(V3r 2 ). Analogously we ob- 
tained 



and 



P=-V3 



V'(r) 



The stress dependent shear modulus is given by 



G'= Cll ~ Cl2 -P: 



(A.5) 



(A.6) 



f(n,) +3 ^). 



(A.7) 



Then, the condition of vanishing stress dependent shear 
modulus coincides with condition J3J in the static limit. 
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